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ABSTRACT 

The  Naval  Postgraduate  School  (NPS)  Lidar  System  was  modified  to  allow 
comparison  of  lidar  profiles  with  radiosonde  profiles  preparatory  to  evaluation  of  the 
lidar  as  a  monitor  for  changes  in  the  boundary  layer.  The  detector  package  was 
modified  to  permit  day  and  night  operation  and  to  reduce  the  overlap-limited 
minimum  measurement  range.  Redesign  of  the  transmitter  beam  expansion  optics 
raised  the  power  threshold  for  internal  optics  damage  while  improving  the  filling  of 
the  18"  output  mirror  and  maintaining  eye  safe  operation.  Tests  of  the  system  on 
completion  of  modification  demonstrated  clear-air  returns  to  ranges  of  1  to  1.5 
kilometers.  Direct  comparison  of  clear-air  lidar  returns  with  radiosonde  balloon 
launches  can  now  be  undertaken  to  evaluate  correlation  of  changes. 
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I.  INTRODUCTION 


A.  BACKGROUND 

A  lidar  system  was  designed  and  built  at  Naval  Postgraduate  School,  Monterey 
California,  to  provide  range  information  to  atmospheric  features  such  as  clouds.  It 
was  originally  intended  to  be  integrated  into  the  Naval  Postgraduate  School’s 
Infrared  Search  and  Target  Designation  system  (NPS-IRSTD).  The  lidar  system 
consisted  of  a  laser  transmitter,  employing  a  frequency-doubled  NdrYAG  laser,  a 
telescope  to  expand  the  beam  in  order  to  meet  laser  safety  requirements,  and  a 
receiver  employing  a  photomultiplier  tube  [Ref,  1].  For  this  thesis  it  was  initially 
intended  to  use  the  Naval  Postgraduate  School  (NPS)  lidar  to  make  a  comparison 
with  radiosonde  balloon  atmosphere  profiles  to  determine  if  the  lidar  system  could 
be  used  to  monitor  changes  in  the  atmospheric  vertical  profile.  The  system  could 
then  be  used  as  a  cuing  device  to  indicate  when  a  new  radiosonde  launch  is  required 
to  maintain  an  accurate  radar  coverage  prediction. 

The  lidar  system  was  noted  to  have  the  following  problems  that  were  required 
to  be  corrected  prior  to  its  use: 

1.  The  detector  was  limited  to  night  time  operation  due  to  high  background  noise. 

2.  The  laser  transmitter  was  burning  the  mirrors  used  to  direct  the  beam  to  the  tele¬ 
scope. 

3.  The  trigger  signal  was  inconsistent. 

4.  The  receiver  sensitivity  was  not  adequate  for  clear  air  lidar  signals. 

5.  The  field  of  view  was  too  small  for  near  field  returns. 

6.  The  detector  noise  was  too  high  to  detect  weak  signals. 

7.  A  ground  loop  problem  between  the  transmitter  and  the  receiver  existed. 

8.  A  power  supply  was  needed  to  supply  plus  and  minus  five  volts,  plus  and  minus 
six  volts,  plus  twelve  volts,  and  minus  eighteen  volts  to  the  detector. 
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B.  RESEARCH  OBJECTIVES 

Due  to  the  number  of  problems  with  the  lidar  system  the  objective  of  this  thesis 
was  changed  from  comparison  of  the  lidar  profiles  to  radiosonde  balloon  data  and 
for  its  use  as  a  cuing  device  to  the  design  modification  of  the  Naval  Postgraduate 
School  Lidar  system  to  enable  it  to  see  clear  air  lidar  signals  out  to  a  range  of  1-2 
kilometers  during  daytime  operation  and  to  eliminate  the  noted  problems. 

After  the  system  is  modified  and  at  a  later  date  it  is  intended  to  perform  the 
comparison  to  balloons  and  test  the  feasibility  of  a  lidar  system  to  be  used  as  a  cuing 
device  for  launches.  The  system’s  initial  design  requirements  to  be  able  to  range  to 
clouds  will  be  mmntained  in  the  modified  system. 

The  system  will  be  required  to  meet  the  laser  safety  requirements  for  eye 
safety.  Use  of  a  system  that  is  not  eye  safe  is  not  allowed  for  obvious  reasons. 


C  LIDAR  SYSTEM  DESIGN  OVERVIEW 

The  modifications  to  the  existing  lidar  system  were  limited  to  those  needed  to 
make  the  system  operational  for  the  purpose  of  measuring  clear  air  returns.  Other 
applications  may  be  desired  and  added  at  some  time  in  the  future. 

The  overall  system  will  remain  basically  the  same  as  it  will  use  the  frequency- 
doubled  NdiYAG  laser  in  the  transmitter  and  the  beam  will  be  expanded  to  an  18 
inch  diameter  for  eye  safety  purposes.  The  laser  safety  calculations  are  detailed  in 
Chapter  III.  The  receiver  consists  of  a  telescope  to  collect  the  backscattered  laser 
signal  and  associated  optics  to  filter  background  noise  and  focus  the  signal  onto  a 
photomultiplier  tube.  The  receiver  design  is  covered  in  detail  in  Chapter  III.  The 
data  collection  was  performed  using  a  digital  oscilloscope  and  stored  on  a  floppy 
disk  for  off-line  analysis.  The  oscilloscope  is  accurate  enough  to  allow  the 
determination  as  to  whether  or  not  a  lidar  signal  is  present. 
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This  thesis  is  organized  as  follows: 


1.  Chapter  II  is  an  Introduction  to  Lidar  Theory  and  the  atmospheric  parameters 
that  affect  radar  propagation. 

2.  Chapter  III  is  the  Design  of  the  Lidar  System  including  laser  safety  require¬ 
ments. 

3.  Chapter  IV  is  Data  Collection  and  Analysis. 

4.  Chapter  V  is  Conclusions  and  Recommendations. 
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II.  LIDAR  THEORY 


A.  INTRODUCTION 

The  use  of  lidar  systems  for  monitoring  the  atmosphere  is  not  a  new  idea.  It  has  been 
looked  at  in  great  detail  and  should  someday  be  able  to  provide  the  same  data  that  currently 
is  only  available  from  radiosonde  balloons.  The  two  primary  methods  being  considered  for 
determining  the  water  vapor  and  temperature  profiles  are;  [Ref.  2] 

1.  Raman  scattering,  and 

2.  differentia]  absorption. 

The  basic  lidar  system  is  made  up  of  three  sections.  The  first  is  the  transmitter  section 
which  includes  the  laser  and  the  associated  optics  to  direct  the  laser  energy  in  the  desired 
direction.  The  second  is  the  receiver  section  which  includes  the  light  gathering  optics, 
detector  and  amplifiers.  The  third  is  the  data  processing  section  which  includes  an 
oscilloscope  and  a  computer  for  signal  analysis. 

This  thesis  is  concerned  with  the  use  of  an  aerosol  lidar  as  a  cuing  device  to  indicate 
when  a  radiosonde  launch  is  needed;  the  theory  of  how  basic  lidar  systems  work  will  be 
covered  in  detail  in  the  following  sections. 


R  LIDAR  EQUATION 

The  amount  of  energy  arriving  at  the  detector  is  dependent  upon  the  initial  laser 
energy  transmitted,  the  probability  that  the  light  will  be  reflected  back  toward  the  detector, 
and  the  probability  that  the  reflected  light  will  strike  the  detector.  The  basic  lidar  equation 
that  describes  these  factors  is  given  as;  [Ref.  3] 

R 

,  ^  -2rK(R)dR 

A  (cx)  J  2-1 


where 


Pj  is  the  energy  received  by  the  detector; 

Pq  is  the  initial  energy  transmitted  by  the  laser; 

Ao  is  the  effective  aperture  of  the  receiver  optics; 

R  is  the  range  from  the  detector  defined  as  cAt/2; 

^(X.)  is  the  transmission  efficiency  of  the  receiver  optics  at  wavelength  X; 

^R)  is  the  geometric  form  factor  ; 

^(X,R)  is  the  backscatter  coefficient  at  R  and  wavelength  X; 

c  is  the  speed  of  light; 

X  is  the  laser  pulse  width; 

k(X)  is  the  attenuation  coefficient  of  the  atmosphere  for  wavelength  X. 

Since  the  wavelength  of  the  laser  and  the  detector  system  are  the  same  for  a 
monochromatic  aerosol  lidar  the  wavelength  dependance  will  be  omitted  from  the 
lidar  equation  for  the  rest  of  this  report. 

The  transmission  efficiency  factor  |(X)  is  a  constant  value  which  is  only 
dependant  on  the  lenses  and  filters  used  in  the  receiver.  The  geometric  form  factor 
KR)  is  comprised  of  two  factors.  The  first  is  the  overlap  function  which  is  broken 
up  into  three  categories  for  a  biaxial  lidar  system.  The  second  is  the  effect  due  to 
the  reflected  lidar  signal  from  the  near  field  not  being  in  the  focal  plane  of  the 
receiver  telescope.  This  results  in  the  received  signal  increasing  as  range  is 
increased  and  the  image  plane  moves  to  the  focal  plane  for  an  object  at  infinity.  The 
signal  has  a  1/R^  dependence  but  because  of  the  geometric  effect  previously 
discussed  the  actual  lidar  signal  increases  until  the  range  reaches  a  value  at  which 
most  of  the  irradiance  of  the  backscattered  lidar  signal  is  focused  onto  the  detector. 
At  this  point  the  lidar  signal  will  decrease  as  1/R^.  The  three  categories  of  the 
overlap  function  are  (a)  no  overlap,  (b)  complete  overlap,  and  (c)  partial  overiap  (see 
Figure  1).  The  first  two  can  be  easily  handled  but  the  third  is  a  difficult  problem. 
[Ref.4] 

Both  components  of  the  geometric  form  factor  are  dependent  on  the  system 
configuration  or  geometry  and  do  not  vary  with  each  laser  firing.  Since  a  change  in 
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a)  NoOveriap 


Figure  1.  Three  overlap  situations  possible  for  a  biaxial  lidar.  R  is  the  radius  of  the 
laser  illumination,  W  the  radius  of  the  receiver  field  of  view,  d  is  the  axis  spacing. 
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the  atmospheric  conditions  is  the  factor  of  concern,  the  geometric  factors  cancel  out 
when  comparing  the  change  from  one  profile  to  the  next  This  allows  the  use  of  the 
partial  overlap  area  to  indicate  a  change  in  atmospheric  factors  even  though  the 
actual  parameters  themselves  would  be  difficult  to  measure.  The  geometry  of  a 
basic  lidar  system  would  depend  on  the  laser  location  with  respect  to  the  detector. 
This  would  include  the  separation  of  the  receiver  and  laser  axes  (d),  the  laser  beam 
divergence  (0)  and  the  receiver  telescope  opening  angle  (())).  Figure  2  shows  the 
basic  geometry  of  a  biaxial  lidar  system.  [Ref.  3] 


Figure  2.  Geometry  of  a  Biaxial  Lidar  System.  W  is  the  radius 
of  the  laser  circle  of  illumination  and  r^  is  the  radius  of 
the  circular  field  of  view. 

The  next  term  in  the  lidar  equation  is  the  atmospheric  volume  backscatter 
coefficient  ^(X,R).  P(X,R)  is  the  backscattered  intensity  per  unit  solid  angle  per  unit 
path  length  (kilometers'^  steradian'^).  P(X,R)  can  be  expressed  as 


P(X,R) 


471 


2-2 


where  N(R)  is  the  number  density  of  the  species  at  range  R  responsible  for  the 
backscattered  radiation  at  wavelength  X  and  a^(X)  is  the  relaUve  cross  section.  Since 
(/(X)  is  a  constant,  the  variation  in  P(X,R)  is  due  to  the  number  density  change.  For 


maritime  areas  the  primary  aerosol  is  water  vapor  which  viiries  as  a  function  of 
altitude.  It  is  this  parameter  that  can  be  used  to  indicate  a  change  in  the  relative 
humidity.  And  again  since  it  is  only  a  change  that  is  to  be  detected  to  indicate 
whether  or  not  a  new  radiosonde  launch  is  needed,  the  change  in  the  lidar  signal 
would  be  an  indication  of  a  change  in  P(X,R),  which  would  be  an  indication  of  a 
change  in  the  relative  humidity. 

The  next  term  in  the  lidar  equation  is  the  transmission  factor  which  can  be 
expressed  as: 

T(KR)  -  -2rK(KR)dR 

J^o  2-3 


Where  k(X,R)  is  the  extinction  coefficient  as  a  function  of  wavelength  and  range. 
The  wavelength  is  constant,  but  the  extinction  coefficient  varies  due  to  variation  in 
the  number  of  absorbers  and  scattering  molecules  as  the  range  changes.  The 
extinction  coefficient  can  be  further  reduced  to 


k(X.,/?)  -  +B(R,X,0)  2A 

Where 

</(X,R)  is  the  atmospheric  absorption  coefficient,  and 
B(X,R,0)  is  the  atmospheric  scattering  coefficient  from  which  P(X,,R)  is 
derived. 

The  wavelength  dependence  in  equation  2-4  will  be  omitted  due  to  the  use  of 

a  mcmochromatic  laser.  The  o^R)  term  is  the  sum  of  the  aerosol  and  molecular  (Mie 
and  Rayleigh)  scattering  coefficients.  The  dominant  absorbers  and  scatterers 
depend  chi  the  frequency  of  the  light  propagating  through  the  atmosphere.  The  range 
dependence  is  due  to  the  variation  in  the  number  density  of  the  absorbers  and 
scatterers  along  the  path  of  the  laser  beam. 
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C  RECEIVER  THEORY 

In  order  for  the  lidar  system  to  provide  useful  information  the  receiver  section 
must  be  able  to  detect  the  backscattered  laser  signal  against  the  background  noise. 
This  section  will  describe  photomultiplier  tube,  pre-amplifier,  and  logarithmic 
amplifier  operation. 

The  photomultiplier  tube  is  a  vacuum  tube  photoemissive  detector  consisting 
of  photocathode,  some  number  of  dynodes,  and  an  anode  as  shown  in  Figure  3. 
When  photons  of  sufficient  energy  strike  the  cathode  free  electrons  are  generated. 
The  electrons  are  accelerated  towards  the  anode  by  an  electric  field  and  dynodes  are 
placed  such  that  as  the  electrons  move  toward  the  anode  they  strike  the  first  of  a 
chain  of  dynodes.  The  dynodes  are  maintmned  at  a  high  potential  and  for  each 
electron  that  strikes  the  dynode  6  secondary  electrons  are  generated.  These 
secondary  electrons  then  strike  the  next  dynode  in  the  chain,  generating  more 
electrons.  The  total  gain  is  dependent  on  the  number  of  dynodes,  the  material  used 

and  the  potential  between  dynodes.  The  gain  can  be  described  as  G=6^  where  N  is 
the  total  number  of  dynodes.  [Ref.  5] 

Some  disadvantages  of  using  photomultiplier  tubes  are: 

1.  High  voltage  is  required  to  achieve  the  desired  gain. 

2.  The  wavelengths  are  restricted  to  the  visible  and  ultraviolet  spectrum. 

3.  The  tubes  are  fragile. 

4.  The  quantum  efficiencies  are  generally  low  (10-25%). 

The  advantages  of  the  photomultiplier  tube  are: 

1.  High  gains  which  allow  detection  of  very  weak  signals. 

2.  Low  internal  noise. 

3.  High  signal  gain  to  noise  gain  ratio, 

4.  Large  detector  area, 

5.  Fast  response  times,  typically  on  the  order  of  2  nanoseconds. 
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Vaoium  Tube 


Figure  3.  Photomultiplier  tube  block  diagram  [Ref.  5]. 


There  are  four  types  of  noise  that  must  be  considered  when  describing  detector 
systems.  These  are: 

1.  Shot  noise  (i{sj).  This  noise  is  the  result  of  the  statistical  generation  of  electrons 
by  the  cathode.  The  rms  value  of  the  shot  noise  current  is  given  as  [Ref.6]: 

'at  -  2-5 

where 

ig  is  the  current  on  the  cathode  from  the  lidar  backscattered  signal, 
i^  is  the  dark  current  (the  thermally  generated  current  in  the  absence  of 
any  light), 

it,  is  the  current  generated  by  the  background  light, 
e  is  the  electron  charge, 

Af  is  the  system  bandwidth, 

G  is  the  gain  of  the  photomultiplier. 

2.  Johnson  or  thermal  noise  (ij).  This  is  the  noise  generated  from  the  thermal  ex¬ 
citation  in  a  resistor.  The  rms  value  of  the  Johnson  noise  can  be  written  as 
[Refs.  6  and  7]: 

hrms  y  /j  2-6 

where  R  is  the  equivalent  resistance  of  the  photomultiplier  and  load. 

3.  Background  ncHse  (i|^.  This  is  the  noise  component  which  results  from  the  ra¬ 
diant  spectral  energy  from  the  sun.  Since  this  noise  is  the  limiting  source  of 
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noise  during  daytime  operation  of  photomuitiplier  tube  detectors  it  will  be  fully 
discussed  here.  The  background  noise  level  can  be  reduced  by  matching  the 
field  of  view  of  the  detector  to  that  of  the  transmitted  laser  field  of  view  as 
closely  as  possibly  and  by  the  use  of  a  narrow  band  filter.  The  background  cur¬ 
rent  can  be  written  as  [Refs.  6  and  7]: 


^  Xer){\)E^(X) 


2-7 


where 

X  is  the  wavelength  of  the  background  radiation, 

it(X)  is  the  quantum  efficiency  of  the  photomultiplier  tube, 

Eb(X)  is  the  background  spectral  energy, 
e  is  the  electron  charge, 
c  is  the  speed  of  light, 
h  is  Plank’s  constant,  and 

Tj  is  the  detection  interval,  which  is  not  the  same  as  the  laser  pulse  width. 

Assuming  that  the  terms  that  describe  the  background  energy  (E^CX),  which  is 
dependent  on  the  wavelength)  are  constant  over  the  bandwidth  of  the  optical 
filter,  Eb(X)^  can  be  written  as: 


where 

AX  is  the  narrowband  filter  wavelength  interval, 

Sb(X)  is  the  background  spectral  radiance, 

to(X)  is  the  transmittance  efficiency  of  the  receiver  optics, 

Qq  is  the  receiver  solid  angle, 

Aq  is  the  receiver’s  effective  aperture,  and 
is  the  detection  interval. 

Substituting  equation  2-8  into  equation  2-7  ^elds  the  following: 

^  XeT)(X)AX5^(X)i^(X)QoAo 
h  "  2-9 

4.  Quantum  noise.  This  noise  is  generated  from  the  statistical  fluctuations  in  the 
generation  of  the  secondaiy  electrons  at  the  dynodes  in  the  photomultiplier 
tube.  When  the  gain  of  the  first  dynode  in  the  chain  is  high,  the  dynode  noise 
is  negligible  and  can  be  neglected.  [Ref.  6] 
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The  signal-to-noise  ratio  (SNR)  can  be  written,  from  equations  2-5  and  2-6, 
as:  [Refs.  6  and  7] 


SNR  - 


4kTAf, 


2-10 


When  R  is  large  enough  the  Johnson  noise  is  small  compared  to  the  shot  noise 
term  and  can  be  neglected.  The  two  cases  left  to  consider  are  daytime  and  night 
time.  These  will  be  discussed  in  detail  since  they  determine  the  maximum  range  of 
the  lidar  system. 


1.  Daytime  Lidar  Range  Determination 

When  the  lidar  will  be  operated  during  daytime  the  background  noise  is  the 
limiting  factor  and  ib  » i*  +  id-  Therefore  the  SNR  equation  can  be  simplified  to: 


SNR  - 


leT^AfG^ 


2-11 


For  single  pulses  the  minimum  detectable  signal  can  be  found  by  letting 
SNR=1.  This  is  the  best  case  situation  and  is  the  theoretical  minimum.  The 
minimum  detectable  power  can  be  determined  by: 

1.  Noting  that  the  signal  current  can  be  written  as  a  function  of  the  incident  power 
P(r):  [Ref.  6] 

.  2.12 

^  he 

2.  Now  substituting  equations  2-7  and  2-12  into  equation  2-11  and  solving  for 
P(r)  with  SNR  =  1  yields: 
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2-13 


P{r) 


2hcAKSf^(k)x^(k)Q^AQAf 


By  solving  the  lidar  equation  (2-1)  for  range  (r)  and  substituting  equation 
2-13  for  P(r)  the  maximum  theoretical  daytime  range  can  be  written  as: 


max 


PqX^  ( r)  JcAjaii^exp  [  -2fo~o(r')dr-] 
JShAkSi,(k)xQ(X)QQAf 


2-14 


where 

fmax  is  maximum  range  of  the  lidar  system  (meters), 

Pq  is  the  transmitted  power  (watts), 

X  is  the  laser  pulse  width  (seconds), 

P(r)  is  the  atmospheric  volume  backscattering  coefficient  at  range  r 

(meters’^  *  steradians'^), 
c  is  the  speed  of  light  (meters  per  second), 

Ao  is  the  effective  aperture  of  the  receiver  (meters^), 

X  is  the  wavelength  of  the  background  radiation  which  due  to  the  narrow 
band  filter  can  be  assumed  to  be  the  same  as  the  laser  wavelength 
(meters), 

Ti(X)  is  the  quantum  efficiency  of  the  photomultiplier  tube, 
a(r)  is  the  atmospheric  extinction  coefficient  (meters'^), 

h  is  Plank’s  constant  (watts  *  seconds^), 

AX  is  the  narrow  band  filter  width  (microns), 

Sb(X)  is  the  spectral  radieuice  of  the  background  (watts  *  meters'^  * 

steradians'*  *  microns"^), 

Xq(X)  is  the  receiver  optics  transmission  efficiency  at  wavelength  X, 

Qq  is  the  receiver  solid  angle  (steradians),  and 
Af  is  the  receiver  bandwidth  (Hertz). 


2.  Night  Time  Lidar  Range  Determination 

During  night  time  operation  the  signal  current  i^  is  much  greater  than  the 

background  and  dark  currents  which  can  therefore  be  neglected.  The  SNR  for  night 
time  can  then  be  written  as: 
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SNR  - 


2-15 


/2g2 

leT^AfG^ 

which  simplifies  to; 

SNR - .  2-16 

2ei^Af 

If  the  photomultiplier  rise  time  is  shorter  than  the  laser  pulsewidth  the  signal 
current  i^  is  approximately  the  same  as  the  average  i^  value  and  equation  2-16  can 
be  further  simplified  to: 


or 


2-17 


The  night  time  maximum  range  can  be  found  by  substituting  equation  2-19  into 
the  basic  lidar  equation  (2-1).  The  receiver  optics  transmission  efficiency  is 
accounted  for  by  multiplying  P(r)  by  to(X).  The  resulting  expression  for  the 
maximum  night  time  range  is  given  by: 
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fmat  is  the  maximum  range  of  the  lidar  system  (meters), 

Pois  the  transmitted  power  (watts), 

T  is  the  laser  pulse  width  (seconds), 

P(r)  is  the  atmospheric  volume  backscattering  coefficient  at  range  r 
(meters  *  *  steradians'*), 

Ao  is  the  effective  aperture  of  the  receiver  (meters  ), 

X  is  the  wavelength  of  the  background  radiation  which  due  to  the  narrow 
band  filter  can  be  assumed  to  be  the  same  as  the  laser  wavelength 
(meters), 

T](X)  is  the  quantum  efficiency  of  the  photomultiplier  tube, 
o(r)  is  the  atmospheric  extinction  coefficient  (meters  *), 
h  is  Plank’s  constant  (watts  *  seconds^), 

Tq(X)  is  the  receiver  optics  transmission  efficiency  at  wavelength  X,  and 
M  is  the  receiver  bandwidth  (Hertz). 

The  primary  reason  for  deriving  the  maximum  daytime  and  night  time  range 
equations  is  to  provide  a  means  to  estimate  the  capability  of  the  system  to  meet  our 
objectives  for  different  atmospheric  conditions. 


D.  RADAR  PROPAGATION  IN  THE  ATMOSPHERE 

Wii  ;n  attempting  to  accurately  model  the  propagation  of  radar  (microwaves)  it 
is  required  to  know  the  radio  refractivity.  The  radio  refractivity  (Nj^j-)  is  described 
by  the  following  equation:  [Ref.  8] 

-  T7.6j  -  5.6  j  +  3.73  x  lO^p  2-21 

where  P[mb]  is  the  total  atmospheric  pressure,  T[K]  is  the  temperature,  and  e[mb] 
is  the  water  vapor  partial  pressure. 

Ducts  are  caused  by  large  changes  in  humidity  and  temperature  which  are 
typically  found  in  the  sea-surface  boundary  layer.  It  is  these  ducts  and  how  they 
have  changed  since  the  previous  radiosonde  launch  that  are  of  particular  interest 

Currently  radiosonde  balloons  are  used  to  measure  the  atmospheric  parameters 
of  interest  in  the  calculation  of  the  refractivity  index.  There  is  however  no  method 
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to  verify  that  the  radiosonde  data  is  still  valid  after  a  few  hours.  Since  a  change  in 
the  primary  parameter  that  is  me.  cured  by  the  radiosonde  (relative  humidity)  would 
also  result  in  a  change  in  the  lidar  return  signal  [Ref.  9],  the  lidar  system  could  be 
used  to  monitor  the  atmosphere  and  indicate  when  there  is  a  need  for  a  new 
radiosonde  launch. 

For  a  lidar  system  to  be  used  as  discussed  above  it  needs  to  be  able  to  detect 
changes  in  the  boundary  layer.  This  requires  that  the  system  be  able  to  look 
vertically  approximately  1-1.5  kilometers. 

There  is  a  strong  interest  in  the  atmospheric  marine  boundary  layer  since  this 
is  where  the  atmospheric  effect  on  radar  and  electro-optical  signal  propagation  is  the 
greatest  It  would  be  a  significant  tactical  advantage  to  know  when  the  atmospheric 
conditions  have  changed  enough  to  significantly  degrade  the  current  performance 
predictions.  This  Thesis  will  concentrate  on  the  modification  of  the  NFS  Lidar 
System  to  allow  clear  air  lidar  returns  out  to  a  range  of  1-1.5  kilometers.  Once  the 
system  has  been  modified  the  comparison  of  the  lidar  vertical  profiles  to  the 
radiosonde  balloon  profiles  can  be  performed. 
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III.  LIDAR  SYSTEM  DESIGN 


A.  OVERALL  DESIGN 

The  initial  design  was  done  to  achieve  the  maximum  possible  range  within  the 
constraints  of  eye  safety  and  the  availability  of  the  laser.  The  initial  lidar  system 
consisted  of  a  frequency  doubled  Nd:YAG  medical  laser,  a  17.75  inch  telescope  to 
expand  the  beam  to  meet  eye  safety  requirements,  a  detector  containing  a 
photomultiplier  tube  and  associated  optics  to  focus  the  signal  onto  the  tube,  a  digital 
oscilloscope  to  collect  the  data  and  a  computer  for  data  storage  and  off-line  analysis. 
A  block  diagram  of  the  initial  design  is  shown  in  Figure  4.  Since  the  system  is  to  be 
used  to  monitor  the  atmosphere  near  the  Naval  Postgraduate  School  grounds  and  the 
airspace  of  the  Monterey  airport  the  main  constraint  is  eye  safety.  [Ref.  1] 

The  basic  overall  design  will  remain  basically  the  same;  however  the  internals 
of  both  the  receiver  and  laser  transmitters  are  significantly  modified  to  allow  clear 
lidar  returns  suitable  for  the  end  use  of  comparison  to  radiosonde  data.  The  specific 
design  changes  will  be  discussed  including  the  reasons  for  each. 


R  LIDAR  TRANSMITTER 
1.  Laser 

The  laser  used  in  the  transmitter  is  still  the  frequency-doubled  Nd:YAG 
laser  model  MK  100  manufactured  by  Kigre  Incorporated.  The  laser  head  and 
power  supply  are  separate  units.  The  power  supply  unit  contains  the  control  section 
for  the  laser.  The  number  of  pulses  per  shot  is  selected  in  steps  of  one  from  1  to  5 
and  the  laser  output  power  has  eight  different  selectable  settings  (see  Table  1).  A 
power  meter  is  also  part  of  the  power  supply  unit.  Table  1  lists  the  laser 
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Laser  and  associated  optics 


Tektronics 
DSA  602A 
Digital 
Oscilloscope 


photomultiplier  mbe 


Hoppy  Disc 
Data  Transfer 


Figure  4.  Initial  Lidar  System  Block  diagram  [Ref.  1]. 
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Table  1:  LASER  SPECIFICATIONS  [REF.l] 


Characteristic 

Value 

Comment 

Maximum  Energy  Output 

mean  - 

9.126  millijoules 

From  Table  1  of  Ref.  1 

(5  pulses  per  shot) 

a  =  1.587  millijoules 

Beam  Diameter 

4.0  +0.2  millimeters 

From  Table  1  of  Ref.  1 

Half  Angle  Beam 

0.165 +.0.001  milliradians 

From  Table  1  of  Ref.  1 

Divergence 

Energies  Available 

Laser 

Measured 

Measured  using  an 

(Single  shot  only) 

Setting 

Values 

RJ7100  power  meter 
with  an  RJ736  probe. 

.9 

.47 

Note:  measurements  were 

1.2 

.68 

taken  with  the  side  of  the 

1.8 

.65 

laser  box  removed. 

2.5 

.912 

4.0 

.87 

The  output  energy  indi- 

5.0 

1.22 

cates  that  the  laser  should 

7.5 

1.36 

be  returned  to  Kigre  Incor- 

10.0 

2.06 

porated  for  repair. 

All  values  in  millijoules 

2.  Dall'Kirkham  Telescope 

The  telescope  was  inspected  and  it  was  decided  that  it  should  still  meet  the 
needs  for  this  thesis  but  should  either  be  replaced  or  refurbished  in  the  future  due  to 
the  pitting  that  was  previously  noted  [Ref.l].  The  Dall-Kirkham  is  a  Cassegrain 
telescope  and  its  specifications  are  listed  in  Table  2. 

The  telescope  has  a  3.125  inch  diagonal  mirror  to  redirect  the  focal  plane 
90  degrees  back  onto  the  telescope  axis.  The  telescope  is  mounted  on  a  Pelco 
FT2000L  Pan  and  Tilt  Unit  to  allow  lidar  ranging  and  profiling  in  azimuth  and  ele¬ 
vation.  The  pan  and  tilt  unit  can  change  elevation  to  +  90^  and  azimuth  0  to  355^. 
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Table  2:  DALL-KIRKHAM  TELESCOPE  SPECIFICATIONS  [Ref.  1] 


Manufacture 

Tinsley  Laboratories  Incorporated 

3900  Lakeside  Avenue 

Richmond,  California  94806 

Primary  Mirror 

i 

Diameter:  17.75  inches  (aspheric) 

Radius  of  Curvature:  126  inches 

Focal  length:  63  inches 

Optical  coating;  Al/SiO 

Substrate;  Cervit 

F  number:  3.55 

Secondary  Mirror 

Diameter  5.75  inches  (spherical) 

Radius  of  Curvature:  44.8  inches 

Optical  coating:  Al/SiO 

Substrate:  Cervit 

Overall  System 

Effective  focal  length;  159.75  inches 

Mirror  separation:  46.25  inches 
(vertex  to  vertex) 

Effective  F  number  9.0 

Focal  plane  is  21  inches  behind  the  vertex 
of  the  primary  mirror. 

3.  Laser  Optics 

The  initial  transmitter  box  was  configured  such  that  the  laser  output  first 
went  through  a  beam  splitter  that  directed  the  1.06  pm  portion  of  the  beam  to  a  beam 

stop.  The  beam  was  then  redirected  90°  by  a  mirror,  and  then  the  second  mirror 
redirected  the  beam  another  90°.  After  the  beam  was  redirected  by  the  second 
mirror  it  was  expanded  by  a  -74mm  focal  length  lens.  The  expanded  beam  was  then 
redirected  out  of  the  transmitter  box  onto  the  3. 125  inch  diagonal  in  the  telescope  by 
a  mirror  attached  to  a  precision  mount  (see  Figure  S). 

This  design  resulted  in  the  first  mirror  being  burned.  This  problem  had  to 
be  corrected  before  the  iidar  could  be  used.  The  first  attempt  to  correct  the  problem 
was  to  install  99.9%  reflective  mirrors  in  place  of  the  less  reflective  ones  initially 


20 


Figure  5.  Initial  Laser  Transmitter  Box  Design, 
used.  This  did  correct  the  burning  problem  but  then  it  was  noted  that  the  beam  was 
expanding  too  much  before  exiting  the  laser  optics  box.  This  was  resulting  in 
approximately  40-50%  of  the  laser  energy  being  lost  with  a  significant  reduction  in 
the  backscattered  laser  signal.  The  other  problem  noted  with  the  transmitter  was  that 
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the  trigger  signal  was  inconsistent.  This  needed  to  be  corrected  before  data  could  be 
collected  with  the  oscilloscope. 

The  laser  beam  was  over-filling  the  secondary  mirror  which  was  further 
reducing  the  return  signal.  Another  problem  discovered  was  that  the  laser  box  was 
about  1/4"  too  far  toward  the  rear  of  the  telescope.  This  resulted  in  the  output  beam 
not  being  uniform. 

Since  several  problems  existed  with  the  transmitter  the  decision  was  made 
to  change  the  configuration  to  correct  the  noted  problems  and  to  provide  easier 
alignment  of  the  transmitter. 

The  first  modification  was  to  expand  the  laser  beam  before  the  first  mirror. 
This  completely  eliminated  the  possibility  of  burning  any  mirrors.  Second  by  using 
a  telescope  to  control  the  beam  expansion  it  was  significantly  easier  to  set  the  beam 
divergence  to  just  fill  the  secondary  mirror. 

Next  the  45°  532  nanometer  reflective  mirror  was  used  to  redirect  the 
beam  through  the  second  lens  of  the  collimating  lens  set  (telescope),  see  Figure  6, 
and  onto  the  3.125"  diagonal  minor.  This  filtered  out  the  1.06  pm  wavelength  and 
reflected  99.9%  of  the  532  nm  wavelength. 

The  PIN  diode  was  moved  out  of  the  laser  box.  A  fiber  optic  cable  was 
installed  to  detect  the  laser  pulse  off  the  rear  of  the  laser  head  and  to  provide  the 
signal  to  the  PIN  diode  which  was  located  in  a  box  attached  to  the  end  plate  of  the 
laser  box,  see  Figure  7.  This  configuration  resulted  in  a  trigger  that  was  not  only 
stable  but  much  larger  in  strength,  ~l-2volts,  than  the  original  design.  The  laser 
trigger  signal  ensured  that  the  data  collection  was  triggered  on  an  actual  laser  pulse 
vice  the  laser  flashlamp  noise. 
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lOX  Beam  Expander 


4.5  inches 


When  used  as  a  lOx  beam  expander 


Figure  6.  Collimating  lens  set. 

The  beam  expanding  optics  were  placed  to  fill  the  secondary  mirror.  The 
radius  of  the  secondary  mirror  was  chosen  as  the  99%  point.  The  spacing  between 
the  first  and  second  lenses  of  the  collimator  was  calculated  in  order  to  provide  the 
required  beam  expansion  taking  into  account  the  distance  to  the  focal  point.  The 
matrix  used  to  calculate  the  solution  for  two  lenses  is: 


3-1 


whe.e 

w’  is  the  beam  radius  on  the  secondary  mirror  at  the  1/e  point  (48.68  millime¬ 
ters). 
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0’  is  the  half  angle  beam  divergence  required  to  fill  the  secondary  at  the  1/e 
point  and  knowing  the  focal  plane  location  (27.53  milliradians), 

T2  is  the  distance  from  the  second  lens  of  the  collimating  set  to  the  secondary 
mirror  (millimeters), 

f2  is  the  focal  length  of  the  second  lens  of  the  collimating  set  (100  millime¬ 
ters), 

Ti  is  the  distance  from  the  first  to  the  second  lens  of  the  collimating  set  (mil¬ 
limeters), 

fj  is  the  focal  length  of  the  first  lens  of  the  collimating  set  (-12.7  millimeters), 
Wq  is  the  laser  beam  waist  radius  (2  millimeters), 

60  is  the  laser  beam  half  angle  divergence  (0.1 65  milliradians). 

The  45°  532  nanometer-reflective  dichroic  mirror  was  placed  0.5  inches 
from  the  first  lens  of  the  collimating  set.  This  prevents  the  beam  from  over  filling 
the  mirror.  The  beam  diameter  is  8  mm  (0.315  inches)  at  the  mirror  which  is  1.25  x 
1.75  inches.  Also  the  3.125  inch  diagonal  is  of  sufficient  size  to  redirect  the  entire 
beam  to  the  secondary  mirror  since  the  beam  diameter  at  the  diagonal  is  38.78  mm 
(1.53  inches). 

The  collimating  lens  set  was  adjusted  to  fill  the  secondary  mirror  by 
mcmitoring  the  beam  at  the  secondary  mirror.  The  distance  from  the  first  lens  of  the 
collimating  lens  set  to  the  secondary  mirror  is  151.76  cm.  The  final  half-angle  beam 
divergence  out  of  the  Dall-Kiikham  telescope  was  measured  and  found  to  be  0.362 
milliradians.  This  would  meet  the  laser  safety  requirements  since  the  laser  beam 
would  never  focus  in  the  far  field. 


C  LASER  SAFETY 

Laser  safety  is  a  major  concern  in  lidar  design  and  operation.  The  laser  beam 
had  to  be  greatly  expanded  to  meet  the  requirements  for  the  maximum  allowed 
energy  density,  the  Maximum  Permissible  Exposure  (MPE),  of  0.5  x  10"^  joules  per 
centimeter^  [Ref.  10].  It  is  assumed  for  laser  safety  purposes  that  99%  of  the  energy 
that  enters  the  telescope  reaches  the  secondary  mirror.  This  should  be  a  reasonable 
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Figure  7.  Laser  Transmitter  Box. 

assumption  since  the  beam  has  been  expanded  to  fill  the  secondary  mirror  (see 
Equation  3-1).  The  secondary  blocks  a  portion  of  the  beam  that  exits  the  telescope 
and  the  reduced  energy  will  be  accounted  for  in  the  laser  safety  calculations. 

The  range  at  which  the  laser  beam  can  safely  be  viewed  is  called  the  Nominal 
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Ocular  Hazard  Distance  (NOHD)  [Ref.  10].  The  geometry  of  the  beam,  the  laser 
power,  and  MPE  all  define  the  NOHD.  The  primary  means  to  prevent  damage  to 
personnel  is  to  ensure  they  are  either  outside  the  NOHD  or  wear  protective  goggles. 
Table  3  lists  the  NOHD  for  all  the  laser  settings.  NOHD  is  calculated  using  the 
equation: 

^  _  Qexp[-\ir]  ^  , 

r(a+r<j))i2 


where 

H  is  the  radiant  exposure  (joules  per  centimeter'^), 

Q  is  the  pulsed  laser  radiant  output  (joules  and  defined  at  the  1/e  point), 
p  is  the  atmospheric  attenuation  co^icient, 
r  is  the  range  (centimeters), 

a  is  the  laser  beam  diameter  at  the  1/e  point  (centimeters), 

<|>  is  the  laser  beam  divergence  angle  (milliradians). 

The  reduction  in  the  output  energy  will  be  corrected  for  in  the  final  calculation 

of  the  NOHD  for  each  laser  setting.  Since  the  laser  beam  is  Gaussian  the  fraction  of 

energy  blocked  by  the  secondary  mirror  can  be  determined  by  [Ref.  11]: 

•—  -  l-exp\ — =— I  -  0.4045  3-3 

Pt  \  dI) 


where 

Pg  is  the  power  blocked  by  the  secondary  mirror, 

Pj  is  the  total  power  of  the  laser  beam, 

Ds  is  the  diameter  of  the  secondary  mirror  (centimeters), 

D^  is  the  diameter  of  the  laser  beam  at  the  1/e  point  (centimeters). 

The  denominator  of  equation  3-2  must  be  corrected  for  the  area  blocked  by  the 

secondary  mirror  also.  The  denominator  can  be  rewritten  as: 
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(Dj^+r^) 

2 


3-4 


Taking  into  account  that  Q  in  equation  3-2  is  the  energy  contained  in  the  beam 
between  the  1/e  points  Q  can  be  written  as  0.86  X  Qp  where  Qp  is  the  total  energy 

in  the  laser  beam  [Ref.  1 1].  Equation  3-2  can  now  be  written  as: 


H 


( 1  -  0.4045)  0.86gy.ejcp[-nr] 
r(Oz. +'•<!>)  *1^  [Dg-f 

2  J  "12. 


3-5 


where 

Dl  =  28.7  cm, 

H  =  MPE  =  0.5  X  10-^  j-cm-2 

Dg  =  14.61  cm 

<J>  =  0.724  X  lO'^rad. 

The  question  of  what  h£q)pens  to  the  energy  tliat  was  blocked  by  the  secondary 
mirror  is  one  of  concern  when  considering  eye  safety.  Some  is  reflected  back  into 
the  transmitter  box  and  some  is  absorbed  by  the  walls  of  the  inside  of  the  telescope 
which  are  painted  black.  The  amount  that  leaves  the  telescope  is  a  matter  that  needs 
further  review.  This  could  be  done  using  a  power  meter  and  comparing  the  output 
energy  with  that  predicted  by  the  laser  safety  calculations.  The  potential  eye  hazard 
can  then  be  determined. 
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Table  3:  NOHD 


Laser  Setting 

Laser  Energy  Measured  at 
the  Laser  Box  Exit 

Qj-  (millijoules) 

NOHD  (meters) 
Assuming  no  atmospheric 
attenuation 

10 

2.06 

347.4 

7.5 

1.36 

219.2 

5.0 

1.22 

190.3 

4.0 

0.87 

110.7 

2.5 

0.912 

120.9 

1.8 

0.65 

53.5 

1.2 

0.68 

61.7 

.9 

0.47 

0.6 

The  worst  case  NOHD  is  calculated  in  Table  3  by  assuming  there  is  no 
atmospheric  attenuation.  The  actual  NOHD  should  be  less  than  the  values  in 
cdumn  3  of  Table  3.  The  laser  output  energy  for  each  setting  is  less  than  the  values 
for  the  same  setting  in  Reference  1.  This  indicates  that  the  laser  needs  to  be  sent 
back  to  the  manufacture  for  repairs.  It  was  also  noted  during  the  firing  of  the  laser 
that  the  power  would  drop  and  become  inconsistent  if  the  laser  was  fired  with  less 
than  approximately  two  seconds  between  shots.  This  should  not  result  in  a  laser 
safety  problem  since  the  energy  went  down. 


D.  LIDAR  RECEIVER 

The  original  receiver  was  not  usable  during  daytime  due  to  noise.  Its 
sensitivity  had  to  be  improved  to  allow  the  monitoring  of  the  boundary  layer,  which 
is  the  final  objective  of  this  thesis.  A  ground  loop  problem  existed  between  the 
transmitter  and  the  detector  which  would  result  in  signal  degradation  or  loss. 
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It  was  decided  that  the  detector  system  should  be  modified  to  correct  the  noted 
problems.  The  first  step  was  to  replace  the  lens  with  a  telescope  to  gather  the 
backscattered  light,  see  Figure  8.  A  suitable  cassegrain  telescope  was  found  with  a 
focal  length  of  750  mm  and  an  F  number  of  6.  The  telescope  was  mounted  on  top 
of  the  receiver  box  that  houses  the  photomultiplier  tube  (PMT).  To  limit  any  stray 
light  from  entering  the  system  the  telescope  was  coupled  to  the  PMT  by  first 

redirecting  the  532  nm  wavelength  signal  90°  with  a  532  nm  45°  selective  mirror. 
This  would  provide  the  first  filter  to  limit  the  background  light.  The  532  nm  light  is 
then  collimated  with  an  F-3.5  lens  to  allow  the  signal  to  be  passed  at  normal 
incidence  through  two  10  nm  wide  filters  centered  at  532  nm.  The  light  shield  tube 

and  housing  for  the  45°  minor  were  painted  black  to  absorb  as  much  stray  light  as 
possible  and  reduce  the  background  noise.  The  use  of  a  tube  prevents  light  leaks  in 
the  box  from  introducing  more  background  noise.  The  F-3.5  lens  was  placed  to  also 
focus  the  signal  onto  the  PMT.  To  provide  the  required  fine  adjustment  capability 

needed  to  ensure  the  signal  strength  is  maximized,  the  45°  mirror  assembly  was 
made  with  vernier  control  of  the  left/right  and  up/down  adjustment  This  was  for  the 
final  alignment  adjustments  after  the  initial  coarse  alignment. 

The  field  of  view  was  not  restricted  with  a  pinhole  since  the  objective  is  to  look 
at  the  boundary  layer.  This  will  result  in  a  reduced  maximum  range  to  targets  such 
as  clouds  for  daytime  operation  only.  The  field  of  view  is  0.5067  x  10'^  steradians 
or  a  half  angle  subtense  of  12.7  milliradians.  The  range  at  which  beam  overlap 
begins  is  found  using  Equation  3-6.  The  range  where  complete  overlap  occurs  can 
be  found  using  Equation  3-7. 


{R  +  l)^  +  RQ  -  0.0889 

3-6 

(/?+  1)<|>  -  /?e  + 0.5389 

3-7 

Where 
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R  is  the  range  (meters), 

(j>  is  the  receiver  half  angle  field  of  view  (milliradians), 

6  is  the  laser  beam  half  angle  divergence  (milliradians). 

Solving  Equations  3-6  and  3-7,  the  initial  beam  overlap  starts  to  occur  at  a 
range  of  5.9  meters  and  complete  overlap  at  42.6  meters  (see  Figure  9). 

The  original  PMT  was  determined  to  be  excessively  noisy  and  was  replaced 
with  a  new  R1913.  An  R928  was  considered  as  a  replacement  but  during  testing 
proved  to  be  less  sensitive  than  the  new  R1913.  The  R1913  PMT  specifications  are 
listed  in  Table  4.  A  LeCroy  Model  VV  1(X)B  wideband  pulse  amplifier  was  installed 
in  the  PMT  box  and  was  coupled  to  the  PMT  to  minimize  the  capacitance.  This  was 
done  to  reduce  the  RC  time  constant  as  much  as  possible  and  at  the  same  time  to  use 
a  500  ohm  resistor  to  prevent  over-driving  the  PMT.  The  amplifier  has  a  gain  of  10, 
a  rise  time  of  less  than  2  nanoseconds,  and  a  bandwidth  of  greater  than  200  MHz 
[Ref.  12].  The  measured  c^)acitance  was  25  pF.  This  provided  a  resolution  of  12.5 
ns  or  1.875  m.  This  should  be  sufficient  for  the  atmospheric  monitoring. 

The  output  of  the  LeCroy  amplifier  was  AC  coupled  to  an  Analog  Modules, 
Model  382  Logarithmic  Amplifier.  The  logarithmic  amplifier  has  an  input  dynamic 
range  of  -200p.V  to  -2  V,  3  nanosecond  risetime,  and  an  output  of  -.5  to  +.6  volts 
[Ref.  13].  The  logarithmic  amplifier  would  compress  the  dynamic  range  of  the 
return  signal  since  it  would  only  amplify  the  weaker  signals.  As  the  range  increases 
the  signal  strength  decreases  and  the  amplification  by  the  logarithmic  amplifier 
increases. 

The  ground  loop  problem  was  corrected  by  making  insulating  sleeves  for  the 
bolts  that  fasten  the  receiver  assembly  to  the  mounting  bracket  and  by  the  use  of 
insulating  spacers  and  washers. 
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Table  4:  R1913  PHOTOMULTIPLIER  SPECIHCATIONS  [Refs.  1  and  14] 


Cathode  Sensitivity 

65  milHAmperes/Watt  (420  nanometers) 

55  milliAmperes/Watt  (532  nanometers) 

Current  Amplification  (Gain) 

3.0  X  10^ 

Rise  Time 

1.0  nanosecond  (2000  VDC) 

Anode  Dark  Current 

Typical  0.05  nanoAmperes 

Maximum  0.5  nanoAmperes 

Spectral  Response 

Range  185-900  nanometers 

Peak  Wavelength  420  nanometers 

Quantum  Efficiency 

Average  in  the  Visible  Region  13% 

At  532  nanometers  15% 

The  use  of  a  PMT,  pre-amplifier,  and  a  logarithmic  amplifier  required  several 
different  voltages  to  supplied.  The  voltages  required  were  +5  V,  +  6V,  +12  V,  and 
-18  V.  Since  a  power  supply  that  could  supply  all  the  required  voltages  does  not 
exist  commercially,  one  was  designed  and  built  (see  Figure  10). 

E.  DATA  ACQUISITION 

The  data  collection  was  initially  performed  using  the  Hewlett  Packard 
HP541  ID  Digitizing  Oscilloscope  which  has  a  sampling  rate  of  Igigasample  per 
second.  Hard  copies  of  the  data  were  taken  using  a  dot  matrix  printer  and  then 
scanning  the  displays  for  inclusion  in  this  thesis.  Additional  data  was  collected 
using  the  Tektronix  DSA  602A.  The  primary  reason  for  shifting  to  the  Tektronix 
was  the  increased  ease  of  data  storage  since  the  DSA  602A  has  a  built  in  disk  drive. 
The  DSA  602A  also  has  a  2  gigasample  per  second  sampling  ability  which  could 
provide  a  better  range  resolution  of  the  display.  [Ref.  15] 

The  next  chapter  will  present  actual  data  taken  and  a  discussion  of  the 
maximum  clear  air  lidar  return  and  hard  target  ranges. 
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IV.  DATA  COLLECTION  AND  ANALYSIS 


A.  SYSTEM  DEMONSTRATION 

The  modified  lidar  system  was  tested  on  2  August  1994  at  1430  hours  to  ensure 
that  it  was  able  to  see  clear  air  lidar  signals.  The  system  was  pointed  at  a  distant 
cloud  bank  and  Figure  1 1  shows  a  strong  lidar  return  out  to  a  distance  of  750  meters. 
This  was  before  the  system  alignment  was  completed,  with  the  High  Voltage  set  at 
1000  Volts  and  no  logarithmic  amplifier.  The  system  demonstrated  that 
atmospheric  profiles  can  be  monitored  with  the  modified  system.  The  weather  was 
clear  with  a  cloud  bank  moving  in  from  the  seaward  side  of  Monterey.  The  lidar 
system  was  operated  by  the  Physics  Department  staff  initially  and  then  the  author 
took  data  with  the  Tektronix  DSA  602A  Digitizing  Signal  Analyzer  described  in 
Chapter  III.  The  author  had  completed  the  required  laser  training  to  allow  the 
operation  of  the  system.  The  primary  purpose  of  the  test  was  to  verify  that  the 
system  could  now  be  used  to  obtain  lidar  profiles  of  the  atmosphere.  The  ultimate 
test  is  to  determine  if  the  lidar  system  is  capable  of  good  lidar  returns  of  the 
boundary  layer.  To  test  the  system  several  returns  were  taken  and  a  selection  of 
these  are  displayed  in  Figures  12  -  14. 


a  RANGE  LIMITS 

There  are  two  range  limits  that  apply  to  the  lidar  system.  The  first  is  the 
maximum  range  for  atmospheric  monitoring.  The  second  the  maximum  range  to 
hard  targets  such  as  small  hemispherical  reflectors.  The  first  limit  is  divided  into 
daytime  and  night  time  maximum  ranges. 

The  maximum  theoretical  daytime  atmospheric  lidar  range  can  be  determined 
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from  Equation  2-14,  using  typical  values  for  the  atmospheric  conditions  and  system 
parameters.  The  theoretical  maximum  daytime  range  is: 


max 


jP^r^  (r)  JcAqXt]  (X)  exp  [-2j^"'‘‘o(r')  dr'] 
^/»AXS^().)To(X)QoA/ 


1 


4-1 


where 

^max  is  maximum  range  of  the  lidar  system  (meters), 

Pgis  the  transmitted  power  (0.4  x  10^  watts,  2millijoules/5  nanoseconds), 

T  is  the  laser  pulse  width  (5  nanoseconds), 

P(r)  is  the  atmospheric  volume  backscattering  coefficient  at  range  r 
(the  values  vary  between  1  x  10'^  and  1  x  lOr^  meters'^  *  steradians'^ 
[Ref.  16].  A  value  of  1  x  10'^  will  be  used  [Ref.  16]), 
c  is  the  speed  of  light  (3  x  10^  meters  per  second), 

Ao  is  the  effective  aperture  of  the  receiver  (0.1 13  meters^), 

X  is  the  wavelength  of  the  background  radiation  which  due  to  the  narrow 
band  filter  can  be  assumed  to  be  the  same  as  the  laser  wavelength 

(532  X  10*^  meters), 

il(X)  is  the  quantum  efficiency  of  the  photomultiplier  tube  (15%), 
o(r)  is  the  atmospheric  extinction  coefficient  (the  values  vary  between 
8  X  10^  and  2  X  lO  ^  [Ref.  17].  A  value  of  5  x  10^  will  be  used), 
h  is  Plank’s  constant  (6.626176  x  10*^  watts  *  seconds^), 

AX  is  the  narrow  band  filter  width  (0.01  microns), 

Sb(X)  is  the  spectral  radiance  of  the  background  (20  watts  *  meters'^  * 
steradians'^  ♦  microns"^  [Ref.  6]), 

T^X)  is  the  receiver  optics  transmission  efficiency  at  wavelength  X(0.25 
due  to  the  use  of  two  532  nm  filters  in  series), 

Qq  is  the  receiver  solid  angle  (0.506  x  10'^  steradians" ,  and 
Af  is  the  receiver  bandwidth  (1  x  10^  Hertz). 

The  resulting  equation  is: 

^max  ”  947.76exp[-5x  10"^rJ  (meters)  4-2 

and  must  be  solved  numerically.  The  resulting  maximum  range  is  676  meters. 

For  night  time  the  theoretical  maximum  range  is  given  by: 


J  (meters).  44 

Again  the  solution  is  found  numerically  and  is  852  meters. 

The  maximum  theoretical  ranges  will  be  compared  to  the  actual  lidar  ranges 
obtained  for  daytime.  Sufficient  time  was  not  available  to  allow  recording  of  night 
time  data. 

The  first  lidar  return  was  in  good  agreement  with  the  expected  daytime 
maximum  range.  However  this  was  at  a  lower  High  voltage  than  the  optimum. 
When  an  average  of  several  shots  is  performed  the  range  increases  as  expected. 

The  time  scale  in  Figure  1 1  is  5  pseconds  per  centimeter.  The  return  was  the 
first  strong  clear  lidar  return  and  the  system  had  not  been  aligned  to  fully  maximize 


Figure  1 1.  Clear  air  lidar  return 
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Figure  12.  Clear  Lidar  Return  After  System  Alignment  Completed. 


the  signal.  After  the  system  alignment  was  ccnnpleted  dear  returns  were  seen  out  to 
between  1200  and  1300  meters  (see  Figures  12  through  14). 

The  range  can  be  increased  by  increasing  the  High  Voltage  to  1200  or  1600 
volts.  The  SNR  can  be  enhanced  by  averaging  over  several  shots  since  the 
improvement  is  proportional  to  the  square  root  of  the  number  of  samples  [Ref.  6]. 
Fcm*  example  16  shots  results  in  an  increase  in  the  theoretical  maximum  range  during 
daytime  to  approximately  2600  meters  and  for  night  time  to  approximately  8200 
meters. 

For  hard  target  ranges  the  target  reflectance  and  cross  sectional  areas  must  be 
known  or  a  reasonable  estimate  made.  The  range  equation  for  hard  targets  is  given 
as:  [Ref.  18] 
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where 

K  is  the  beam  profile  function, 

Pq  is  the  transmitted  power  (watts), 

r  is  the  target  laser  cross  section  (meters^), 
is  the  receiver  aperture  diameter  (meters), 
r]x  is  the  transmitter  optical  efficiency, 
r]j  is  the  receiver  optical  efficiency, 

9  is  the  laser  beam  width  (radians),  and 
Pj.  is  the  received  signal  power  (watts). 

The  beam  profile  for  Gaussian  beam  is  given  by: 

A:(tp,  qp)  -  2exp(~^'' ^  ]  4-6 

where 

qp  is  the  laser  beam  width  (radians), 

'll!  is  the  line  of  sight  pointing  error  (radians), 
r  is  the  range  to  the  target,  and 
(0  is  the  Gaussian  waist  radius. 

The  target  laser  cross  section  for  a  sphere  is  given  by:  [Ref.  18] 

r  -  Jipz^  4-7 


where 

r  is  the  target  laser  cross  section  (metera^), 
p  is  the  reflectance  of  the  sphere’s  surface,  and 
z  is  the  radius  of  the  sphere  (meters). 

For  the  NPS  lidar  system  an  estimate  of  a  hard  target  maximum  range  is  made 
using  the  following  assumptions: 

1.  The  hemispherical  reflectance  is  0.191  [Ref.  18]. 

2.  The  targets  radius  is  0.254  meters. 

3.  The  pointing  error  is  1  milliradian  (assumed  value). 
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4.  The  receiver  optical  efficiency  is  25%. 

5.  The  beam  pro^le  function  value  is  assumed  to  be  4,  which  is  the  value  for  a 
uniform  beam  and  is  more  conservative  than  for  a  Gaussian  beam  [Ref.  18]. 

6.  The  transmitter  optical  efficiency  is  60%. 

Solving  for  the  daytime  maximum  range  yields  4700  meters.  If  the  laser  power 
were  to  be  increased  to  9  millijoules  the  maximum  daytime  range  would  be  6850 
meters  which  is  still  less  than  the  required  15  kilometers  for  integration  into  the 
IRST  system.  To  meet  this  requirement  would  require  an  output  power  of  320 
millijoules.  This  would  result  in  a  NHOD  of  8530  meters.  This  could  be  a 
significant  problem  to  friendly  forces. 

Possible  improvements  are; 

1.  Be  able  to  change  the  field  of  view  of  the  receiver  depending  on  the  range  of 
interest.  Hard  target  and  atmospheric  monitoring  are  the  two  cases  of  concern. 

2  Remove  the  second  filter  when  using  the  reduced  field  of  view. 
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Figure  13.  Lidar  and  Hard  Target  Returns. 
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Figure  14.  Lidar  return  and  the  edge  of  a  cloud  bank  rolling  in  from  sea. 
The  next  chapter  discusses  the  conclusions  and  recommendations  for  future 
work- 
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V.  CONCLUSIONS  AND  RECOMENDATIONS 


The  current  performance  level  of  the  telescope  falls  somewhat  short  of  the 
requirement  for  a  small,  unattended,  monitor  of  atmospheric  profile  stability.  Trade 
offs  may  be  made  in  the  design,  but  the  eye  safety  and  range  requirements 
necessitate  a  large  optics  transmitting  apertuie,  dominating  the  size  of  the  system. 
The  NFS  Lidar  System  was  modified  to  enable  the  monitoring  of  the  atmospheric 
marine  boundary  layer.  The  improvements  in  the  system  are: 

1.  The  detector  is  now  able  to  operate  during  daytime  and  provide  the  needed 
range  to  allow  the  boundary  layer  monitoring. 

2.  The  shorter  range  for  beam  overiap  allows  the  monitoring  of  the  boundary  lay¬ 
er  from  approximately  5.5  meters  to  1.5  kilometers. 

3.  The  installation  of  a  beam  expander  in  the  transmitter  has  eliminated  the  prob¬ 
lem  of  burning  mirrors  and  made  it  easier  to  prevent  the  over  filling  of  the  mir¬ 
rors  in  the  transmitter  telescope. 

4.  The  alignment  of  the  laser  beam  has  been  made  significantly  easier. 

5.  The  ground  loop  problem  was  eliminated. 

6.  The  trigger  signal  is  much  larger  and  very  reliable. 

There  are  more  improvements  that  could  be  made  in  the  future.  Some  of  these 
are: 

1.  The  mounting  of  the  detector  to  the  telescope  needs  to  be  modified  to  allow  a 
more  precise  left/right  and  up/down  adjustment. 

2.  A  larger  aperture  telescope  could  be  used  on  the  receiver  system  to  increase  the 
ranges. 

3.  A  new  transmitter  telescope  could  be  used  to  improve  the  transmitter  optical 
efficiency. 
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4.  A  method  of  attenuating  the  energy  that  is  blocked  by  the  secondary  mirror 
needs  to  be  developed. 

5.  The  pan  and  tilt  head  assembly  needs  to  be  modified  to  allow  the  mounting 
bolts  to  clear  the  top  of  the  cover  when  in  the  full  vertical  position. 

Overall  the  system  appears  to  be  capable  of  use  to  monitor  changes  in  the 

atmospheric  marine  boundary  layer.  As  can  be  seen  from  the  figures  in  Chapter  IV 
the  system  can  achieve  clear  air  lidar  returns  out  to  the  required  distance  of  1 
kilometer.  Once  the  pan  tilt  head  is  modified  an  actual  vertical  profile  can  be 
obtained  and  the  atmospheric  structure  should  be  detectable.  It  should  prove  to  be 
a  good  cuing  device  as  to  when  a  new  radiosonde  balloon  launch  is  indeed 
necessary.  The  system  is  now  capable  of  collecting  the  needed  data  to  perform  the 
comparison  of  radiosonde  balloon  profiles  to  lidar  profiles. 

This  thesis  has  been  a  challenging  and  rewarding  experience  in  lidar  design.  It 
has  given  the  author  a  better  understanding  of  the  problems  associated  with  the 
research  and  design  of  a  system. 
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